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Abstract 



i <~| The 125 GeV resonance discovered at the LHC could be a heavy quarkonium pseudoscalar meson. The 

diagonalization of the mass matrix of the isoscalar quarkonium states \Q and |£° ) produces the states 
identified with the pseudoscalar heavy quarkonium mesons C° and C°'. For a mixing angle <j> ~ 36 ° 
the mass of the bound state pseudoscalar resonance £° is m ? o ~ 125 GeV. The decay rates of the 
quarkonium £° meson are estimated and compared to the standard model Higgs boson predictions. 
The importance of determining the spin-parity of the 125 GeV resonance is discussed. Criteria for 
experimentally discriminating between the pseudoscalar, J PC = f° quarkonium meson and the 
scalar J PC = ++ Higgs boson are investigated. 
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O 1 Introduction 

The search for the Higgs boson at the Large Hadron Collider (LHC) has discovered a new boson. The CMS [1] 
and ATLAS £3] experiments based on 5fb _1 data show an excess of events at ~ 125 — 126 GeV that was 
£N| . evident already in the 2011 data. When combining the 7 TeV and 8 TeV data, both experiments separately 

reach a sensitivity corresponding to a local significance of ~ 5c The resonance could be the discovery of the 
elusive Higgs boson of the Standard Model (SM) [3J. However, we must make certain that the new resonance 
is the Higgs boson that was postulated to generate the electroweak symmetry breaking, which gives masses 
to the SM vector bosons W and Z and to the fermions. The LHC experiments will eventually have the 
sensitivity to accurately determine the decay rates for the final states 77, ZZ* , WW*,Zj, bb, gg and t+t~ 
predicted by the Higgs boson production mechanism. The excess of events at ~ 125 GeV is mainly driven by 
the "golden" decay channels 77 and ZZ* — > 4 leptons. In particular the latter decay channel is a significant 
signature of the Higgs coupling to a vector boson. However, the present sensitivity for this decay channel 
is about 3c There are anomalies in the decay channel observations, such as the zero or negative signal for 
t+t~ decay and that the signal size for the diphoton decay channel is about 1.5-2 times that predicted for 
the Higgs boson. Moreover, the data for the bb decay channel is consistent with the predicted Higgs boson 
rate. 

In the following, we will attempt to identify the new boson with a neutral, heavy quarkonium bound 
state called (° [4]. The standard quantum chromodynamics (QCD) quark-gluon force binds the quarks and 
anti-quarks. An important test of the possible identification of the new boson as the C° is its spin and parity, 
J PC = K Due to the observed decay of the new boson into a diphoton state, it has to have either spin 
or spin 2. Because the SM elementary Higgs boson is necessarily a scalar J PC = ++ boson with even 
parity, it is critically important to determine the parity of the new boson. This can be done at the LHC by 
analyzing the correlations of the angular distributions of the 77 and 4-lepton decay states of the 125 GeV 
boson. 
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2 Mixed Quarkonium Resonance 

The bottomnium and unobserved toponium are isoscalar states \B) = \bb) and \T) = \tt) of heavy quarko- 
nium. With respect to an effective interaction Hamiltonian, heavy quarkonium appears in two different 
isoscalar states |C°) and |C° )• The effective Hamiltonian is given by 

Heft = T~Lq + % mass • (i) 

where 

n mass = K T MK. (2) 

Here, M. is the mass matrix: 

M=( m <' m «' ), (3) 

where |£) and |£') are states of quarkonium that interact through the mixing contributions to^< and K = 
^ I . The mass matrix can be diagonalized: 

D = R T MR, (4) 

where R is the rotation matrix: 



sin m cos < 



and I? is the diagonal matrix: 



Here, bit ~ 2toj ~ 346 GeV and mj ~ 2to;, ~ 9 GeV where we have used the measured quark masses: 
m t = 173 GeV and m b ~ 4.5 GeV. 

In terms of the states \T) and \B), we have 

H' mass = m T TT^ +m B BB\ (7) 

where 

B y V C / V c °s </>C — sin 0C / 
From we obtain the masses: 

my = cos 2 (fim^i — 2 sin <^> cos 4>m^Qi + sin 2 to^ , (9) 

and 

ttlq = cos 2 </>to.£ + 2 sin</>cos ifim^^ + sin 2 to,£<. (10) 

By inverting (fj|, we get 

M = RDR T (11) 

which leads to the result 

_ / m^/ \ _ / cos 2 0tot + sin 2 0tot3 cos <^> sin 0(mB — mj) \ 

~~ \ m cc m C / \ cos</)sin0(TOB — tot) cos 2 0to_b + sin 2 0tot /' 



We derive from (TT2"]) the results 
and 

The off-diagonal term is given by 



TO£ = cos 2 ^>tob + sin 2 0tot, (13) 
m^' = cos 2 (firriT + sin 2 0ms- (14) 
to<^/ = cos0sin0(TO73 — tot). (15) 
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We can determine the mixing angle <f> from the equation: 

4> — arccos[(mT — m c)/ ( m T — ttt-b)] 1 ^ 2 ■ (16) 

The mixing angle cf> ~ 36 is obtained from (JTS]) and from (fT3"|) and (fH)) . we get the masses of the quarkonium 
states \Q and |C') : m (° ~ 125 GeV and m^ a > ~ 230 GeV. We identify the new boson resonance discovered 
at the LHC with the £° bound state quarkonium resonance. The quarks and anti-quarks are bound together 
by the QCD gluon force with a corresponding binding energy Eb- A possible non-relativistic heavy quark 
gluon potential is 

V = --+br, (17) 
r 

and the QCD coupling constant a s (Mz) = 0.118 [5]. 

The short life-time of the top quark for the decay t — > bW + , r t ~5x 10~ 25 sec, results in the toponium 
not forming a bound state. However, the toponium state exists as an state of heavy quarkonium. The heavy 
C° with a mass m, ' ~ 328 GeV does not form a detectable bound state. 



3 Prominent Decays of £° and Higgs Boson 

The (° meson will be primarily produced at the LHC in gluon-gluon fusion with the largest cross section for 
the 1 Sq ground state, and this state will decay into gg,jj,cc,bb, t + t~ , /i+ ^ ,WW* and ZZ*. At leading 
order and in the narrow-width approximation, the production cross section for a boson B in pp collisions is 
given in terms of the gluon decay width by 

2 rl 

<j(pp^B + X) = £^T(B^gg) J dx^-g(x,Q 2 )g(T B /x,Q 2 ), (18) 

where g(x, Q 2 ) is the gluon parton distribution function, and tb = m's/s where s is the pp collision energy 
squared. We identify B with either the pseudoscalar meson B — (° or the Higgs particle B = H°. 

We can obtain approximate results for the £° boson decay widths, using the colored non-relativistic heavy 
quark model [6] . The quark wave function for the lowest lying S-wave state is 

V>(r) = -Liio(r). (19) 

V47T 

In the Coulomb approximation, we have 



1^ ,. ., 

rriQO 



= C 4 a s 3 m c o, (20) 



where Cj is the quark charge and color factor. The parametrization of the bound state wave functions, their 
binding energy and the R{0) at the origin have been computed in ref. [?]■ 
The form for the two-particle £° annihilation decay width is [8j |9] 

T(C° -+ AB) = 327t2(1 \ 6ab) N 1/2 ^ </ml mym^ E W 2 " (2D 



in 



C 



Here, M are the amplitudes and AB — gg, 77, ZZ*, WW*,Z^, bb, cc and 8ab = 1 for the cases gg, 77, ZZ*. 
Moreover, N(x,y,z) — x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. The production cross section in proton-proton 
collisions for £° — > 77 is given by 



a(pp -> C° ^ 77) = ^3 r (C° -► 77) jf dx~g(x, Q 2 )g(r/x, Q 2 



(22) 



where BR((° — > 77) is the branching ratio. The ratio of the C° — > 77 production cross section to that of the 
SM model Higgs boson H is 

u aiC " ± 77) (23) 
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For the mass ttlq = m# ~ 125 GeV the ratio jj, could be greater than unity in agreement with the current 
excess signal for the diphoton decay of the new boson [TJ [2] . 

The M for the W + W~ final state is given in the heavy quarkonium zero velocity approximation v ~ 

by 

M X \C° -> W+W~) = -lw W ~ l ~ M (hoho + (-) X hx)\g 2 cos 2 9 - 2(3 2 vl 2 v S X0 S- X0 g 2 cos 2 8 g^- . (24) 

where A, A are specific helicities of the W bosons, which can take on the values 0,±1. Moreover, 9 is 
the t-channel scattering angle, g and sn\0\y are the SM gauge parameters, = m^/mw and (3% = 

y/l - K/ffl() 2 . 

The M for the ZZ final state is given in the zero velocity approximation v ~ by 



M^(C° ^ ZZ) = -7z _|A| - |X| («ao<Jao + HXa) 1 



cos Ow 



^9 2 [[\-\ sin 2 9 W J cos 2 + - sin 4 9 W 



2g 2 mj 



cos 2 6w 



(cos 2 - § sin 2 9 W ) 2 cos 2 9 sin 2 (9 

o 2 2 9 9 

zmj> — m z — m z 



(25) 



where jz — m^/mz and (iz = yl - (mz/m^) 2 . 
Finally, for the Z7 final state we get 



|M(C° Zi) | 2 = 3 ^^" g2 f 1 cos 2 - I sin 2 6 W V, (26) 



cos 2 6*vy V ^ 3 

where Q g is the quark charge and only A^ = A 7 = ±1 helicity values are permitted. 

As in the SM with a Higgs boson, the coupling of the £° to Z and W bosons is stronger than the coupling 
to photons, and we have that r(C° -> ZZ*) > T{H° -> 77). 

For the S-wave initial £° state and v — > 0, we have for the squared amplitudes for the gg and 77 final 
states: 

]T |M(C° -> gg)\ 2 = S gg Q q (4ira s ) 2 (27) 

spins, color 

and 

J2 |M(C° -> 77)| 2 = ^ 77 (47ra) 2 , (28) 

spins 

where Sj are spin and color factors. We find from the ratio of (l27l) and (|28|) that r(£° — > gg) 3> r(£° — > 77). 
For the decay £° — 77, we get to leading order in the Born approximation: 

r(C° -> 77) = C 77 f — ) |i? (0)| 2 = C 77 a 2 a s 3 m c o. (29) 

We can compare the leading order partial decay width of (° — > 77 to the partial width of the scalar Higgs 
particle decay in the light mass Higgs limit [TU1 EJ H2] : 

where Gf — 1.6637 x 10~ 5 GeV -2 is Fermi's constant. In the limit of a Higgs mass mn ~ 125 GeV and for 
the W and top quark loop contributions, I(w+t) ~ — 10, we obtain 

T(H° 77) - 31keV. (31) 

For the gluon final state, we get for the £° meson decay: 

r(C° gg) = \c gg (-£-) \r (o)\ 2 = \c gg a s 5 m Q0 . (32) 
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The partial decay width in the Born approximation for the decay H° — > gg is given by 

Measurements of the total production cross sections <tq and <jh and the branching ratios of decay channels 
at the LHC can be used to distinguish between the C° pseudoscalar heavy quarkonium resonance and the 
SM Higgs boson. 

4 Parity in the Diphoton Decay of (° 

Let us consider the non-relativistic quarkonium state as the product of the orbital wave function and the 
spin vector S [IB] : 

* n .l, m (r\S,S z ). (34) 
The spin vectors are linear combinations of products of the individual quarks: 

|S = l,S» = l) = |t>|t>, 

|s = i,s, = o) = -^(|t>U> + U)|t>), 

\S = l > S z = -l) = \i)\i), 

|S = > £z = 0) = -^(|t>U>-U)|t», (35) 

where | f)| f) is the product of an anti-quark spin state with S z — +1/2 and a quark state with S z — +1/2, 
producing an state of total spin 5 = 1 and z component S z = I. In analogy with the states of positronium, 
the first three expressions above describe the possible three states of "ortho-quarkonium" , while the last 
expression describes the spin state of "para-quarkonium" . 

By applying the parity operator to the quarkonium orbital wave function, we get 

P*n,l,m{ r ) = *n,i,m(-r) = (-l)'* n ,/, m (r). (36) 

The intrinsic parity of a quark is the opposite of the anti-quark. The overall parity of a quarkonium state 
is —(—1)' and since both ortho- and para-quarkonium decay from a I — state both have odd —1 parity 
at decay. Since parity is conserved in electromagnetic interactions, the photon states must have odd parity. 
The wave function of the photon with momentum k and transverse polarization e is 

A(r,t) = exp(i(k-r- urb)\e). (37) 

The symmetrized wave function for the two-photon state with momenta k and — k and polarizations ei and 
€2 is given by 

A(r,t) =e3q>(i(k.(r-2wi)|ei)|e 2 )+e!xp(i(-k.(r-2wt)|e2)|ei>, (38) 
where r = ri — r 2 and |ei) and |e 2 ) are orthogonal. This expression can be rewritten as 

A(r, t) = -[|ei)|e 2 ) + |e 2 ) | ei )]{exp(i(k • r - 2wi) + exp(i(-k ■ r - 2ojt)} 

+ \[h)\e 2 ) - |e 2 )|e 1 )]{exp(i(k • r - 2wt) - cxp(z(-k • r - 2wt)}. (39) 

In the last equation the terms both exhibit Bose symmetry of the interchange of the Bose indices, 1 < — > 2, 
and the first term has even and the second term odd parity. The para-quarkonium decay state is described 
by the second term in Q: 

A PQ (r,t) = _[| ei )| ea ) - |e a )|e 1 )]{exp(t(k ■ r - 2ut) - eatp(*(-k ■ r - 2wt)}. (40) 

This polarization state is a superposition of different polarization states of the individual photons. Nonethe- 
less, the final measured value of the polarization state is fixed by conservation of parity and the orthogonality 
of the polarization vectors. 
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5 Parity of Decays to ZZ 



The scalar character of the Higgs boson with J PC = ++ can be tested at the LHC. Below the threshold 
for two real Z bosons, the Higgs particle can decay into real and virtual WW* and ZZ* pairs. The partial 
decay width for the Higgs decay into ZZ is given by |15j : 

r(H°^ZZ) = ^^(l-4x + 12x 2 )^/jr^j, (41) 

where x = m 2 z /m 2 H . The parity of the SM Higgs particle is necessarily P = +1. However, a pseudoscalar 
coupling to a Higgs-type particle can be achieved in two-doublet models, in which the coupling to ZZ is 
induced by higher-order loop effects |14j . Moreover, supersymmetry models have multiple Higgs bosons of 
which one can be a pseudoscalar boson A. 

The way to test the parity and spin of the 125 GeV boson is to determine the angular correlations in 
the boson decay to fermion and gauge boson pairs. The observation by the CMS and ATLAS groups of the 
decay or fusion process into two photons rules out the J = 1 directly by the Landau- Yang theorem [16[ 117] 
and determines that C = +1. This would also follow from the decay of the boson into WW* and ZZ* . The 
longitudinal wave function of a vector boson grows with energy in contrast to the transverse wave function. 
The Z boson in the S-wave Higgs production process becomes asymptotically longitudinally polarized. By 
contrast, the Z boson in the £° production process becomes transversally polarized at high energies. This 
experimental signature can be checked at the LHC. The distribution of lepton pairs in the Z — > t~£ + rest 
frame relative to the flight direction of the Z is given by sin 2 9 for longitudinally polarized Z bosons, while 
for transversally polarized states, after averaging over azimuthal angles, it behaves as (1 ± cos 2 9) 2 . 

The azimuthal angular correlations as well as the polar angles depend sensitively on the spin-parity 
assignments of the boson particle. We will study the scalar J PC = ++ Higgs and pseudoscalar J PC = h 
C boson decays [H Hi H3 [21] : 

B^ZZ* ^^£+)(£^£+), (42) 

where B = [H°,£ ]. The reconstructed Z bosons are back-to-back in the B boson rest frame. The two Zs 
are labeled Z\ and Z<i- The reconstructed masses of the Z bosons are denoted by m\ and m%, and they 
define the longitudinal boosts from the B rest frame to the rest frames of the decaying Z\ and Zi bosons: 

and 

7 2 = f^fl-^M). (44) 

We denote by 6±, 4>i and 62, (f>2 the £± and 1% decay angles in the Z\ and Zi rest frames, respectively. 
The decay amplitudes are defined in terms of the two boost parameters 71 and 72: 

7a =7i72(l + /3i/3 2 ) (45) 

7& = 7i72(/3i+/3 2 ), (46) 

where the /8»,7i are the velocities and 7 factors of the on/off-shell Z bosons. The 7 a and 7b satisfy the 
condition 

7a 2 " 7b 2 = 1. (47) 

We have ^ 

la = 5 [ml ~ (ml + ml)}. (48) 

Effective Lagrangians for the HZZ and (ZZ couplings are 

Chzz = {V2G F ) l/2 m 2 z HZ^Z^ (49) 

and 

C C zz = hm 2 z CZ^Z^. (50) 
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Here, Z^ — d^Z v — d v Z^, h is a coupling constant, Z^ v = e fJ - l,pa Z pcr , and ([50")) is a P-wave coupling, odd 
under parity and even under charge conjugation. 

The angular distribution of the leptons in H° — > ZZ — > 41 decay, for on or off-shell Zs of mass mi and 
mi for the scalar Higgs boson is given by |21) : 



dT(H° ->• ZZ*) 2 2 4 h , „„„2/l „„„2/l i I 2 . 5» • 2 



<i cos Videos 02 c 



oc mlm2m H [l + cos 0i cos 02 + (t& + cos <f>) sin 0i sin 



2 , 



72 



+27 a cos sin 0i sin 02 cos0i cos 02 + 2r/ 2 (cos0i cos 02 + 7a cos0sin0i sin 02)]. (51) 
Here, 77 denotes the quantity occurring in the SM couplings of the Z bosons to the final state leptons: 

For the pseudoscalar C° decay C° -> ZZ* -t (£-£+)(£-£+), we obtain [21]: 
dT(C° -> ZZ*) 



d cos Old cos 02C?0 



oc mim27h(l + cos 0i cos 2 - cos ^ sin 0i sin^ 2 + 2rj z cos0i cos0 2 ). (53) 



The general feature that distinguishes the scalar Higgs boson decay from the pseudoscalar £° decay is 
that for scalar particles the production amplitude ~ e% ■ 62 is non-zero for parallel polarization vectors, while 
for the pseudoscalar £° the amplitude ~ ei x e*2 is non-zero for perpendicular polarization vectors. The zero 
spin of the new boson at 125 GeV can be checked experimentally by observing the lack of any correlation 
between the final and initial state particles. 



6 Conclusions 

We have developed a model in which the newly discovered boson at the LHC with a mass 125 — 126 GeV can 
be identified with a heavy quarkonium pseudoscalar resonance C°- By mixing the heavy quarkonium states 
\(°) and |C° ) through a rotation angle 4> ~ 36 °, we obtain two states |C°) and |C° ) with the masses m^o ~ 125 
GeV and m^ ' ~ 230 GeV. The £° can form a bound state with the standard QCD gluon interaction. The 
rates of decays into light hadrons and leptons predicted by the (° can be compared to the decay rates of the 
Higgs boson, so that it can be verified which boson has been observed at the LHC. Moreover, the important 
signature of the differences in the parity of the £° and Higgs boson can eventually be determined at the LHC, 
when enough data and a large enough luminosity have been obtained. This will unambiguously distinguish 
between two of the possible identifications of the new boson as a heavy quarkonium meson or a SM Higgs 
boson. 
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